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We analyze the effects of adiabatic ramping of optical lattices on the temperature of spin-1 bosons 
in a homogeneous lattice. Using mean-field theory, we present the isentropes in the temperature- 
interaction strength (T, Uo) plane for ferromagnetic, antiferromagnetic, and zero spin couplings. 
Following the isentropic lines, temperature changes can be determined during adiabatic loading 
of current experiments. We show that the heating-cooling separatrix lies on the superfluid-Mott 
phase boundary with cooling occuring within the superfluid and heating in the Mott insulator, and 
quantify the effects of spin coupling on the heating rate. We find that the mean-field isentropes for 
low initial entropy terminate at the superfluid-Mott insulator phase boundary. 
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I. INTRODUCTION 

Ultra-cold atoms in optical lattices offer the possibility 
of realizing various fundamental models of strongly cor- 
related bosons and fcrmions A crucial aspect of the 
optical lattice system is its flexibility in controlling differ- 
ent lattice parameters and particle interactions, thereby 
facilitating progress towards the creation of quantum 
emulators. After the observation of the superfluid to 
Mott insulator transition with spin-0 bosons steady 
progress has been made towards trapping spinful atoms. 
These models, along with a precise knowledge and tun- 
ability of the microscopic Hamiltonian can potentially 
lead to a better understanding of quantum magnetism, 
and related phenomena. Unlike magnetic traps which 
freeze the component of spin, optical traps can confine 
^^Na, ^^K, and ^''Rb with hyperfine spin F = 1. Several 
theoretical studies have focused on spinor condensates 
in an optical lattice 0-01 and how the spin degree of 
freedom modifies the phase diagram and the nature of 
superfiuid-Mott insulator transition. Experiments have 
also explored the properties of spinful bosons in harmonic 
traps[8| and, recently, in a double well optical superlat- 
tice 0. 

The technological breakthrough in cooling to ultra-cold 
temperatures paved the way for the realization of Bose- 
Einstein condensation and optical lattice experiments. 
The temperature of a bosonic gas in a trap can be mea- 
sured accurately. However, no established temperature 
measurement exists for optical lattice systems, although 
several proposals have been made [13, El . This makes it 
difficult to obtain a quantitative description of the var- 
ious low temperature phases and thermal and quantum 
phase transitions between them [l2| . 

In current experiments, ultra-cold atoms are first 
loaded in a harmonic trap, and then an external sinu- 
soidal potential created by interfering lasers are slowly 
ramped up to create the optical lattice. For gradual 
enough ramping, there is no heat exchange with the en- 
vironment, and this process can be considered adiabatic 
(constant entropy) [l|, [l^ . It is of great interest to the 



experimentalists to know how the system cools down or 
heats up during the adiabatic process. The change in 
temperature with adiabatic ramp-up of optical lattice for 
the spinless Bosc-Hubbard model and the Fermi Hubbard 
model has been studied by several authors p^ - fisj . 

In this paper, we present the isentropes for spin-1 
bosons in a homogeneous optical lattice for both ferro- 
magnetic and antiferromagnetic spin couplings. We in- 
vestigate the effects of adiabatic ramping on temperature 
within the mean-field approximation, and show that for 
spinor bosons cooling occurs in the superfluid phase and 
heating occurs in the Mott insulator or normal phase. 
We find that the heating-cooling separatrix lies along the 
superfiuid-Mott phase boundary. As the magnitude of 
spin coupling increases, the rate of temperature change 
decreases in the Mott regime, and can both increase and 
decrease in the superfiuid regime depending on the value 
of the spin coupling. We find that the mean-field isen- 
tropes for low initial entropy terminate at the superfluid- 
Mott insulator phase boundary and argue that this is a 
consequence of the absence of breaking of the degeneracy 
of the ground state in the mean field approximation. 

The paper is organized as follows. We discuss the spin- 
1 Bosc-Hubbard model and the details of the mean-field 
theory in section II. In section HI, we investigate the en- 
tropy for this model, and present our results for the isen- 
tropes and temperature changes for different couplings. 
A summary is contained in section IV. 



II. SPIN-1 BOSE-HUBBARD MODEL AND 
MEAN-FIELD THEORY 

The Hamiltonian for spin-1 bosons in an optical lattice 
is given by 

H = -t [ata-j^ + a]^a,<,) + ^Yn^{n^- 1) 

+^E(^'-2-0-E^»"- (1) 
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FIG. 1: (color online). Entropy and filling as functions of 
chemical potential at zero temperature for Uo — 12zt (we 
set zt = 1 in our calculation) and (a) U2 = 0, (b) U2 ~ 
-O.lf/o and (c) U2 = O.lC/o- For U2 = 0, the entropy is 
ln(3) at p = 1, ln(6) at p = 2, and zero elsewhere. This 
can be understood in terms of the degeneracy of spin-1 (three 
component) bosons and total spin at each site. For U2 < 0, 
the energy is minimized with maximum total spin F = 2, and 
hence the entropy at p = 2 changes to ln(5). Similarly, for 
U2 > 0, the ground state has lowest total spin, F — 0, and 
hence the system has zero entropy for p = 2. 



Here a}^ (aia) are boson creation (destruction) opera- 
tors at site i with spin component a {a = 1,0,-1). The 
first term in Eq.([T|) describes the spin dependent hop- 
ping between near neighboring sites. In the second term, 



Uo is the on-site repulsion, and = o-l^o-ia counts 
the total number of bosons on site i. In the third term, 
U2 is the spin dependent interaction which can be zero, 
positive or negative, and Fi = „i al^Faa'dicr' is the 
total spin on site i, where F„„i are the standard spin-1 
matrices. 

The spin dependent interaction, U2, for the spin-1 
model greatly modifies its physics compared to that of 
spin-0 bosons From the symmetry of the bosonic 
wavefunction, scattering with total spin = 1 is prohib- 
ited. The difference in scattering lengths in the F = 
(ao) and F = 2 (02) channels is responsible for the spin 
dependent coupling. The interactions can be expressed as 
Uo = 4:7:h^{ao+2a2)/UI and U2 = 4,T:h^{a2-ao)/3M, M 
being the mass of the atom. The spin dependent interac- 
tion is ferromagnetic when U2 < (02 < ao) and antifcr- 
romagnetic when U2 > (02 > oq). ^"^Na atoms are ferro- 
magnetic and ^'^Rb antiferromagnetic. Our study in this 
paper considers both signs of the interaction. The cou- 
pling constants obey the constraint —1 < U2/U0 < 1/2. 

The zero-temperature phase diagram for the spin one 
model has been calculated with numerical methods such 
as QMC [i,[ll, DMRGpOl, and also within the mean- 
field approximation 0, l2l|. As for the spin-0 case, mean- 
field theory for the spin-1 Bose-Hubbard model 0,113, 113 
captures the essential physics of this system. A finite 
temperature extension of the mean-field analysis was pre- 
sented in 0, revealing the rich phase diagram that in- 
cludes both first and second order transitions. The spinor 
Bose-Hubbard model, with filling of one boson per site 
and for small hopping, can be mapped onto the F = 1 
bilinear biquadratic Heisenberg model, which has been 
studied by many authors |24| - [2q to gain an understand- 
ing of different Mott phases. Here, we will extend the 
mean-field calculation presented in Ref.0 to obtain the 
entropy. 

In the mean-field approximation, the hopping term is 
decoupled as 



(2) 



neglecting fluctuations around the equilibrium value. 
Here, we define ipa = {'^]a) ~ ('^jo-), for a = 1,0,-1 
to be the superfluid order parameter. The use of Eq.([21) 
allows us to rewrite Eq. ([T]) as a sum of independent single 
site Hamiltonians, H = J^i Hf^^ where 



nr = 
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Here we set = 1, where z is the number of nearest 
neighbors. To perform the mean field calculations, we 
write the matrix elements of the Hamiltonian 7?™^ in 
the occupation number basis n^^o, 'T-i.i), and trun- 

cate the onsite Hilbert space Hi by allowing a maximum 



number of particles per site, N„ 
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FIG. 2: (color online). Traversing a first order (a) and second 
order (b) phase boundary in the (/^,T) plane. T — 0.1 (a) 
and T = 0.2 (b). The discontinuity in the entropy in (a) 
satisfies the Clausius-Clapeyron equation relating the entropy 
jump to the slope of the first order phase boundary. Here 
Uo = 12, U2 = 0. 

truncation effects are negligible. We use ^max = 6 for 
simulations with higher number density such as in Fig. [T] 
Since we arc treating a homogeneous system, the site in- 
dex i can be ignored. We diagonalize the Hamiltonian 
to obtain the energy spectrum and eigenstates \4>a), 
and evaluate the partition function and the free energy, 

Z(Ai,C/o,t/2,T;V'.)=5]e-^°/^ (4) 

a 

and 

T{fi, Uo, [/2, T; ^P,) = -Tin Uo, U2, T; V.) (5) 

We set fcs = 1 throughout our calculation. For given 
/i, Uo, U2 and T, the superfluid order parameters ipa are 
obtained by minimizing the free energy, i.e. by solving 
dF/dil}^ = for cr = 1,0,-1. Solving for the extrema 
of the free energy is equivalent to the self-consistency 
condition for ip^. 

After determining the values of the superfluid order pa- 
rameters V'l'i , V'o'^ and ip'^'^ that minimize the free energy, 
other physical quantities can be obtained easily from the 
resulting global minimum of free energy, J-''^ , and eigen- 



values (E^) and eigenstates {4>^) at the global mini- 
mum. The superfluid density is given by, 

and the number density is, 

p = -dTldi, = i ^ e-^»''/^(C''|n|0- (7) 

OL 

Finally, the entropy is calculated from 

S = -dTldT = lnZ+^Yl ^^S^e-^S"/^. (8) 

a 

Pai et al. give a detailed analysis of the superfluid or- 
der parameters ipi^tpo, tp-i for different phases and show 
that in the antiferromagnetic (polar) superfluid, the pos- 
sibilities are ipi = ip-i > 0, V'o = and ipi = = 
Q,ipo > 0. For the ferromagnetic superfluid, ipi — ip^i, 

i>o = V^ipi- 

We use the superfluid density to determine the phase 
diagram in the (T, Uo) plane, ps > corresponds to 
the ferromagnetic or polar superfluid, and ps = cor- 
responds to the Mott insulator phase. At finite temper- 
ature, there is a crossover from the Mott phase to the 
normal phase depending on the value of the compress- 
ibility. 

III. ENTROPY AND ISENTROPIC CURVES 

Figure [TJa) shows the entropy (S) and density (p) at 
zero temperature for U2 = 0,Uo = I2zt (we set the 
energy scale zt ~ 1 throughout our calculation) for in- 
creasing chemical potential along a vertical cross section 
in the (p, Uo) phase diagram. This trajectory traverses 
the superfiuid regions with non-integer p as well as Mott 
plateaux with density fixed at integral values p = I and 
2. We observe that at T 0, although the entropy (S) 
is zero in the superfluid (SF) region, it is nonzero in the 
Mott insulator (MI) lobes. For p = 1, the entropy is 
ln(3), and for p = 2 the entropy S = ln(6). This nonzero 
ground state entropy can be understood as follows. In 
the mean- field treatment, the system Hamiltonian is a 
sum of single site Hamiltonians and in the MI phase for 
one particle per site, three degenerate spin components 
a = -1, 0, 1, contribute to S* = ln(3). For the p = 2 Mott 
phase, the two spin-1 bosons on a single site have total 
spin F = 0, 2, with F ~ 1 eliminated by the symmetry 
constraint on the spin functions. The number of degen- 
erate components is therefore 5 due to F = 2, and 1 due 
to the spin singlet F = 0, so that S = ln(6). 

The entropy for the ferromagnetic case U2 < 0, U2 = 
-O.lt/o, J/o = 12 at T = is shown in Fig. \V[h). A 
negative value of the spin coupling, U2, favors maximal 
total spin. Hence, F — 2 for p = 2, and S = ln(5) is 
reduced from its Uo — value. For p = 1, i.e. the first 
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FIG. 3: (color online), (a) Isentropic curves overlaid on the 
finite temperature phase diagram for C/2 = and density, 
p = 1. During adiabatic ramping of the lattice, the system 
cools along the isentropic lines in the superfluid region, and 
heats in the Mott region. In our mean field analysis, the 
separatrix of cooling and heating lies on the superfluid-Mott 
phase boundary in the (T, Uo) plane, (b) Shows entropy as a 
function of Uo at constant temperatures for intervals of T = 
0.02. The termination of the isentropic lines at the phase 
boundary can be understood by the absence of low entropy 
states as Uo increases into the Mott regime. 



Mott lobe, the entropy is still ln(3) as before, and zero 
in the superfluid regime. Antiferronetic values, U2 > 0, 
favor the F = singlet phase for p = 2 and, therefore, 
5 = for the second Mott lobe. Fig. [IJc) shows the 
entropy for U2 = O.lC/o, Uq = 12 a.t T = 0. As in the 
earlier cases, the entropy is zero in the SF phase, and 
ln(3) in the first MI lobe. The presence of finite entropy 
for the MI phase in mean-field theory even at T = 0, 
influences the topology of the isentropic curves, as we 
shall discuss. 

In contrast to spin-0 bosons, spin-1 bosons in a lattice 
exhibit both first and second order phase transitions at 
finite temperature. Here we examine the entropy across 
these different phase boundaries. In Fig. [2ja), discon- 
tinuities in entropy, density, and superfluid density are 
evident across a first order phase boundary in the (/x, T) 
plane at a constant T = 0.1, and C/q = 12, U2 = 0. Such 
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FIG. 4: (color online). Isentropic curves for a noninteger fill- 
ing, p = 0.7 and U2 = 0. As in the commensurate filling 
case, the heating-cooling separatrix lies on the phase bound- 
ary. However, the low temperature isentropes do not termi- 
nate, and remain under the phase boundary. 



first order finite temperature phase transitions should 
follow the Clausius-Clapcyron equation from thermody- 
namics which relates the discontinuity in the entropy and 
other parameters to the slope of the phase boundary. We 
have verifed that our results indeed satisfy this identity. 
For example, in Fig. [2][a), we find AS/Ap = 3.10 which 
equals the slope dfi/dT in the {fJ.,T) plane. 

As we increase the temperature, the entropy jump de- 
creases, and finally vanishes at the tricritical point. The 
entropy is continous afterwards across a second order 
phase transition as we show in Figl^Jb) for T = 0.2. 

To investigate further the effect of adiabatic ramping 
of the optical lattice, we calculate the entropy at dif- 
ferent couplings (fixed U2 and many Uq) and tempera- 
tures, and construct the isentropic curves in the (T, Uq) 
plane. Fig. EJa) shows the isentropic curves for U2 = 
for fixed occupation number p = 1. Isentropic curves are 
overlaid on the finite temperature phase diagram where 
the boundary separates p^ > (SF phase) and p^ = 
(Ml/Normal phase). During adiabatic ramping of the 
lattice, the system heats or cools following one of these 
isentropic lines. For entropy 1.13 at J/q = the tempera- 
ture starts at T = 0.43 but decreases along the isentrope 
as Uq rises. At the SF boundary, T begins to rise again. 
In fact, all the mean-field isentropes that we obtain show 
this pattern - there is cooling in the SF regime and heat- 
ing in the MI regime with the heating-cooling separatrix 
exactly on the phase boundary. 

The occurence of cooling in the SF phase and heat- 
ing in the Mott phase, and the location of the heating- 
cooling separatrix exactly on the phase boundary can be 
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FIG. 5: (color online). Isentropes overlaid on the finite temperature phase diagram for ferromagnetic (172 < 0) and antifer- 
romagnetic (IJ2 > 0) interactions. The chemical potential jj, is chosen so that density p = 1. For small values of U2, (a) 
U2 = — O.lC/o and (c) U2 = O.lC/o, the isentropes are not much different from each other and from U2 = (Fig. 3). For larger 
values of U2, (b) U2 = —OAUo and (d) U2 = OAUo, the difference is more visible. Along an isentropic line, first we have cooling 
followed by heating for large Uo . The cooling and heating sepatrix again lies on the superfluid-Mott phase boundary. 



understood with the following physical argument - mov- 
ing away from the phase boundary towards higher Uq, 
one enters the Mott phase with reduced number fluc- 
tuations and an integer number of atoms per site. So, 
as Uq increases at constant T, the entropy is reduced. 
Therefore, if we want to keep the entropy constant as Uq 
increases, the temperature must rise, and there is heat- 
ing in the Mott phase. Similarly, moving away from the 
phase boundary towards the SF regime by decreasing Uq, 
more and more particles enter the condensate, reducing 
the quantum depletion. More particles entering the con- 
densate means entropy is decreasing since the condensate 
carries no entropy. So, to keep S constant as Uq decreases 



in the SF, there must again be heating in the system. 

The fact that within mean field theory phase bound- 
aries also demark the switch between cooling and heat- 
ing isentropes is a phenomenon that has also been ob- 
served in studies of classical models of nuclear magnetism 
[13. Recently, it has also been shown that the spin-one 
Blume-Capel model exhibits this same behavior [28j. In 
both cases, it is found that in exact Monte Carlo calcula- 
tions for the same models the heating/cooling separatrix 
continues to track the phase boundary qualitatively, but 
increasingly breaks away as the temperature increases. 

Another feature of the isentropic lines is that for low 
initial entropy (which is also at low temperature) they 
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terminate at the SF-MI phase boundary without ever 
entering the MI phase. Fig. [U showed that within site- 
decoupled mean-field theory, the entropy in the Mott 
regime with p = 1 is ln(3) = 1.0986 even at T = 0. 
As a consequence, if we follow an isentropc with initial 
entropy less than ln(3) in the SF phase, as Uo increases 
beyond the critical value, the isentropc cannot reach the 
MI phase because of its high ground state spin entropy. 
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FIG. 6: (color online). The effect of spin coupling U2 on the 
rate of temperature change. As we vary U2, the slopes of 
the isentropes in the (T, Uo) plane change. Here we show the 
rates of change of temperature with Uo for two isentropes in 
the first Mott lobe and superfluid regime. In the MI region, 
the rate decreases as we move away from U2 = 0; while in the 
SF, it increases or decreases depending on the value oi Uo 

Fig. [Hb) shows entropy as a function of Uq at con- 
stant temperature intervals. As Uo increases into the 
MI at constant temperature, the low entropy values go 
to S* > 1.0986. This illustrates why only the isentropes 
with S > 1.0986 in Fig. [Sja) continue to the MI phase. 
This termination of the isentropic lines is a shortcoming 
of the mean-field theory treatment. There has not been 
any exact numerical study of isentropic lines for the spin- 
1 model within QMC or DMRG. 

In the Mott phase with independent sites, any of the 
three choices ct = 0, ±1 is possible on each site. However, 
it is clear that this degeneracy is broken in perturbation 
theory in t. Two adjacent sites with the same a will have 
a second order lowering /S.E^'^^ = —%t^/U while adjacent 
sites with different a will have a second order lowering 
A£;(2) = -At^/U. The entropy per site is lower than 
ln(3). Thus while mean- field theory successfully captures 
the zero temperature phase diagram as has been verified 
with QMC |6|, and predicts many qualitatively correct 
results for the finite temperature phase diagram, it can- 
not capture the low temperature entropy curves in the 
Mott region. 

Next we examine the isentropes at a non- 
commensurate filling. Fig. m shows the isentropic 
curves for p = 0.7 overlaid on the finite temperature 



phase diagram. As in the commensurate filling case, 
here we also find that the heating-cooling separatrix 
is on the SF-Normal phase boundary. However, the 
low temperature isentropes here do not terminate at 
the phase boundary. Since the phase boundary for 
nonintcgcr fillings never touches the T = line, the 
isentropes can remain under the phase boundary for 
arbitrary C/q- 

We present isentropes for U2 7^ at commensurate fill- 
ing p = 1 in Fig. [S] which shows the isentropic curves for 
both ferromagnetic and antifcrromagnetic cases for two 
different spin couplings. We also show the SF-MI phase 
boundaries on the same plot. The central observations 
are that the isentropes here have the same general prop- 
erty as for f/2 = - there is cooling in the SF phase and 
heating in the MI phase, the phase boundary coincides 
with the heating-cooling separatrix, and the low entropy 
lines do not enter the Mott phase. For U2 = — 0.1[/o and 
U2 = O.IC/qj the isentropes are not much different from 
each other and from the C/2 = case. However for larger 
magnetic coupling, U2 = —QAUq (Fig^Kb)), we find that 
adiabatic ramping results in a slower rate of cooling and 
heating. Starting with the same initial temperature, and 
ending up at the same final optical lattice depth in the 
MI phase, the temperature is lower for U2 = —QAUq than 
for U2 = — O.IJ/q. This is also true for the polar phase 
U2 = 0.4[/o in FiglHd). 

In Fig. [HI we show the slope of the isentropes dT/dUo, 
i.e. the heating rate in the Mott and SF phases, as a 
function of U2. We show results for two isentropes, S = 
1.6 and S = 1.2 in the MI and SF regime. We see that 
the heating rate in the MI regime continually decreases 
as the magnitude of U2 increases. In the SF region, the 
rate is negative denoting cooling. For the two isentropes 
we show, the characteristics of the rates differ based on 
whether we are in a low or higher entropy curve. For 
the species of atoms used in current experiments, ^'^Na 
and ^'^Rb, |J72| is on the order of 0.03. Our mean-field 
study indicates that the effect of magnetic couplings of 
this magnitude on the adiabatic heating and cooling rates 
is small. 

From Fig. [HI we see that the rate of heating in the first 
Mott lobe is maximum at J72 = 0, and slowly decreases 
as magnetic coupling is turned on for both ferromagnetic 
and anti-ferromagnetic case. As noted earlier, increasing 
Uo in the Mott region reduces fluctuations, creating order 
and thereby reducing entropy. Thus the system must 
heat as Uo is increased to keep the entropy constant. Now 
a U2 7^ 0, U2 will be trying to establish another kind of 
order. For example, in the first Mott lobe, U2 > tries 
to establish a bond order in one dimension, and U2 < 0, 
a ferromagnetic order. These other orders compete with 
the simple Mott insulator, slow down the reduction of 
fluctuations as Uo is increased, and result in a slower 
heating rate. 

In the spinless Bosc-Hubbard model, the SF-MI tran- 
sition in the p = 1 and p — 2 Mott lobes are qualita- 
tively similar as they are both a second order transition. 
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FIG. 7: (color online). Isentropes overlaid on the finite tem- 
perature phase diagram for ferromagnetic {U2 < 0) and an- 
tiferromagnetic {U2 > 0) interactions for density p — 2. The 
same values of U2, U2 = — O.lC/o and U2 = O.lUo were used 
as in the p — 1 results in Fig. 4. There are several differ- 
ences from the case of unit filling. Here, the sign of dT/dUo 
can change within the superfluid. The isentropes no longer 
exhibit pure cooling. In addition, the antiferromagnetic isen- 
tropes in (c) have no terminations on the phase boundary. 



The presence of magnetic coupling in the spin-1 model 
changes this scenario For U2 > 0, the SF-MI transi- 
tion in the p = 2 Mott lobe becomes a first order transi- 
tion. Phase boundaries for the even and odd Mott lobes 
also change as a function of U2 d, 0] ■ 

Our final results address the isentropes for density 
p = 2. In Fig. El wc choose U2 = 0, C/2 = -O.lf/o, and 
U2 = Q-lUo- Here also we see cooling in the SF and 




FIG. 8: (color online). The effect of spin coupling U2 on 
the rate of temperature change in the second Mott lobe (p = 
2). Here we show the heating rate for two isentropes S=2.5 
and S=1.85. Unlike the p — 1 case, the heating rate is not 
maximum at U2 = 0, and increases for U2 > 0. 



heating in the Mott regime. Turning on the ferromag- 
netic coupling in Fig. [Tfb) does not have a huge impact 
on the isentropes. However, for antiferromagnetic cou- 
pling U2 = O.lC/o in Fig. [7l[c), we see that the isentropes 
are visibly different as is the shape of the phase bound- 
ary. Inside the SF, the system cools down very slowly. 
Fig. [5] shows the heating rate dT/dUo for a constant en- 
tropy (S=2.5) in the second Mott lobe as U2 is turned 
on. Unlike the case for the first Mott lobe (Fig. [S]), there 
is no maximum at U2 = 0, and the heating rate increases 
monotonically with U2- For U2 > and p = 2, the 
SF/Mott phase boundary is first order, and a spin sin- 
glet phase is formed in the Mott region. The appearance 
of order in the singlet phase greatly affects the heating 
properties. 



IV. CONCLUSION 

In this paper we have extended previous mean field 
theory treatments of the spin-1 Bose Hubbard model to 
compute the isentropic lines. These quantities are of con- 
siderable experimental importance to understand ther- 
mometry in ultra-cold atomic gases in an optical lattice. 

We presented the isentropes in the temperature- 
interaction strength (T, Uq) plane for ferromagnetic, an- 
tiferromagnetic, and zero spin couplings. Following these 
isentropic lines, temperature changes can be determined 
during adiabatic loading of spin-1 bosonic atoms in an 
optical lattice. The isentropes have a number of interest- 
ing features. First, they exhibit pure cooling within the 
SF and heating in the normal/Mott phase. The phase 
boundary precisely corresponds to the location of the 
change in sign of dT/dUo, i.e. the heating-cooling sepa- 
ratrix lies on the SF-MI phase boundary. This can be un- 
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derstood in the symmetric view of heating as one moves 
away in Uq from the phase boundary in either direction 
of SF/Mott. The system gets more ordered with reduced 
entropy. And therefore, temperature must increase to 
keep the entropy constant. Second, the low entropy (low 
temperature) isentropes terminate on the phase bound- 
ary because of the nonzero ground state entropy which 
mean field theory gives for the Mott phase. For noncom- 
mensurate fillings, the isentropes show similar character- 
istics except that the low entropy lines do not terminate 
at the phase boundary. 

The effect of ferromagnetic and antiferromagnetic cou- 
plings on isentropes has been examined, and we quantify 
the rate of heating as U2 changes. We find that in the 
experimentally relevant regime of \U2\ = 0.03 for ^"^Na 
and ^''Rb, the changes in heating-cooling rate is not very 



different from the U2 = case. It would be interesting to 
extend QMC simulations of the magnetic and superfluid 
properties on the spinor Bose-Hubbard model to study 
these thermodynamic properties as well. 
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